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, _ , Abstract 

We have reconsidered theoretical upper bounds on the scalar boson masses within 
■ the two-Higgs-doublet model (THDM), employing the well-known technical condition 

i of tree-level unitarity. Our treatment provides a modest extension and generalization 

of some previous results of other authors. We present a rather detailed discussion of 
the solution of the relevant inequalities and offer some new analytic formulae as well 
as numerical values for the Higgs mass bounds in question. A comparison is made 
with the earlier results on the subject that can be found in the literature. 
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p^; 1 Introduction 

D , The two-Higgs-doublet model (THDM) of electroweak interactions is one of the simplest 
^ ' extensions of the Standard Model (SM). It incorporates two complex scalar doublets in 
the Higgs sector, but otherwise its structure is the same as that of the SM. Obviously, 
^ ■ such a theory is rather appealing on purely aesthetic grounds: in view of the familiar 
. doublet pattern of the elementary fermion spectrum, one can speculate that an analogous 
organizational principle might work for the "scalar Higgs matter" as well. Further, any 
Higgs sector built upon doublets only is known to preserve naturally the famous lowest- 
order electroweak relation p = 1 (where p = m^/(m| cos^ dw))^ which has been tested with 
good accuracy. On the phenomenological side, an important aspect of the THDM is that its 
Higgs sector may provide an additional source of CP violation; in fact, this was the primary 
motivation for introducing such a model in the early literature on spontaneously broken 
gauge theories in particle physics jT]. Of course, there is at least one more reason why the 
THDM has become popular^ during the last two decades or so: its Higgs sector essentially 
coincides with that of the minimal supersymmetric SM (MSSM), but the values of the 
relevant parameters are less restricted. The spectrum of physical Higgs particles within 

^For useful reviews of the subject see e.g. P], [3] 
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THDM consists of five scalar bosons, three of them being electrically neutral (denoted 
usually as h, H and A^) and the other two charged (H^). At present, some partial 
information concerning direct experimental lower bounds for the Higgs masses is available, 
coming mostly from the LEP data (cf. 

On the other hand, it is also interesting to know what could be possible theoretical 
limitations for masses of the so far elusive Higgs particles within such a "quasi-realistic" 
model. For this purpose, some rather general methods have been invented, based mostly 
on the requirements of internal consistency of the quantum field theoretical description of 
the relevant physical quantities. One particular approach, which is perhaps most straight- 
forward in this regard, relies on perturbative unitarity of the S'-matrix. In its simplest form 
it is implemented at the lowest order, by imposing unitarity constraints on the tree-level 
amplitudes of a suitable set of scattering processes. Let us recall that this technique was 
originally developed by B.W. Lee, C. Quigg and H. Thacker (LQT), who employed it in 
their well-known analysis of perturbative upper bound for the SM Higgs boson mass jH]. 
The LQT method was subsequently applied also to electroweak models with extended 
Higgs sectors; some results can be found under refs. 0, [ZI, [HI- In particular, authors 
of the papers [7], jS] analyzed in this way a restricted version of the THDM with CP- 
conserving Higgs sector and obtained slightly differing values of the bounds in question 
(due to slightly different implementations of the LQT method). Recently, the issue of 
tree-unitarity constraints for THDM Higgs boson masses has been taken up again in the 
work jH] (see also Pni,^^)' where a rather general model involving CP violation has been 
considered; this seems to be another vindication of the persisting interest in the subject. 

The purpose of the present paper is to supplement and extend the existing results 
concerning the THDM Higgs mass upper bounds. We carry out a rather detailed analysis 
of a relevant set of inequalities that follow from the requirement of tree-level unitarity. In 
particular, the procedure of explicit solution of these constraints is discussed in considerable 
detail and, among other things, some results of the corresponding numerical calculations 
within a general THDM are presented. For the model without CP violation we were 
able to find a set of analytic expressions as well. Note that in this latter case, most of 
the calculational details are contained also in an earlier unpublished work by one of us 
(see Let us also remark that there is no substantial overlap of the material presented 

in jnmnHH] with our results, so we believe that it makes sense to offer our detailed analysis 
as a contribution to the current literature on the particular problem in question. 

The plan of our paper is as follows: In Sect. |21the THDM scalar potential and the scalar 
fields are described in some detail, in Sect. El we summarize briefly the LQT method and 
its implementation within THDM and in Sect. I^the relevant inequalities expressing the 
tree-unitarity constraints are examined. The main analytic results for the mass bounds in 
question are contained in sections El 13 and Sect. |H1 contains numerical results obtained 
in the CP- violating case (where we have not been able to find analytical results). The 
main results are summarized in Sect. El 



2 



2 THDM scalar potential 



The most general scalar potential within THDM that is invariant under SU(2) x U{1) can 
be written as (cf. [T3j or 0) 



A, 



V{<^) = Ai ($1$! - f 

<l>I<l>i)(<l>^2'^'2)-($l$2)($^3>i) 



A2 ( $5$2 - ^ 



A3($I$l-| + $5$2-|'' 



+A5 



Re($l$2) - ^cos^l +A6 [lm($l$2) - ^sin^ 



Note that such a form involves CP violation, which is due to 7^ [21 • It also possesses 
an approximate discrete Z2 symmetry under $2 —^2', this is broken "softly", by means 
of the quadratic term 



V1V2 (A, 



i.c; COS ^Re(<I>l<I>2) + Ag sin^Im(<l>J$2 



1^1 ^2 Re (Ascos^ - zAesin^) $J$2 (2) 



Let us recall that the main purpose of such an extra partial symmetry within THDM is 
to suppress naturally the flavour- changing processes mediated by neutral scalar exchanges 
that could otherwise arise within the quark Yukawa sector ^3]. Note also that if such a 
symmetry were exact, there would be no CP violation in the Higgs sector of the consid- 
ered model. For further remarks concerning the role of the Z2 symmetry see e.g. jH] and 
references therein. As a quantitative measure of the Z2 violation we introduce a parameter 
z/, defined as 



'XlcosH + Xlsin'^ (3) 

(note that our definition of the u differs slightly from that used in jH].) The minimum of 
the potential (jT} occurs at 



$1 





V2 \vi 



V2 \'^2j 



(4) 



where we have adopted, for convenience, the usual simple choice of phases. Such a min- 
imum determines vector boson masses through the Higgs mechanism; in particular, for 
the charged W boson one gets = ^g^{vf + ?;|), with g standing for SU{2) coupling 
constant. In a standard notation one then writes Vi = vcosj3,V2 = f sin/3, where v is the 
familiar electroweak scale, v = {GpV^)"^^'^ = 246 GeV and /? is a free parameter. THDM 
involves eight independent scalar fields: three of them can be identified with the would-be 
Goldstone bosons , z (the labelling is chosen so as to indicate that they are direct coun- 
terparts of the massive vector bosons W^, Z within an i?-gauge) and the remaining five 
correspond to physical Higgs particles — the charged and the neutral ones h, H, A^. 

We will now describe the above-mentioned Goldstone and Higgs bosons in more detail. 
To this end, let us start with a simple representation of the doublets, namely 



$1 



w 



V2 



1 

hi 



$5 



Wo 



■^{q'^V2 + h2 + iZ2, 



(5) 



(1) 
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Of course, the scalar fields introduced in (0) are in general unphysical; the 2 ^ire taken to 
be complex and the remaining ones real, but otherwise arbitrary. Note that an advantage of 
such a parametrization is that the form of the quartic interactions is then the same as in CP- 
conserving case. The proper Goldstone and Higgs fields are found through a diagonalization 
of the quadratic part of the potential (^. When doing it, a convenient starting point is 
a slightly modified doublet parametrization 



W'r, 



that is obtained from ^ by means of the unitary transformation = /i2Cos^ + 2:2 sin ^, 



= Z2 cos^ — h2 sin^ a w'^ = e ^^wf . Next, the scalar fields in (jH)) are rotated pairwise 



Z2 
as 



H'\ _ / cos/3 sin/3\ fhi\ fA'\ _ / cos/3 sin/3^ ( zi 

h' J y— sin /3 cos (3 J yi^J \ ^ J \~ ^in (3 cos (3 

C\ _ f cos/3 sin/3\ f'wi 

w J sin /3 cos /3 J ytf 2 



(7) 



When the quadratic part of is recast in terms of the new variables, one finds out that 
the z, are massless Goldstone bosons and the represent massive charged scalars. 
At this stage, the fields h', H', A' are still mixed and their mass matrix reads 

I f s2^(Al+A2)+c2g(^c|A5+s|A6) S2/3[-2c2Ai+2s2A2+C2/3(c|A5+s|A6)] |C2;3S2{(A6-A5) \ 

- I S2^[-2c2Ai+2s2A2+C2,3(c|A5+s2A6)] 4 [c^ Ai+s* A2+A3+C2 (c^ As+s^ As)] |s2^S2e(A6-A5) | (9) 

\ |c2/3S25(A6-A5) §S2/3S25(A5-A6) s|A5+c|A6 / 

By diagonalizing it, one gets the true Higgs bosons h, H, A^. The operation of charge 
conjugation C means the complex conjugation of these physical fields (i.e. not of those 
appearing in the parametrization (jSJ). However, we can employ the representation (jH)) 
involving fields that are linear combinations of real variables without complex coefficients. 
Note that for ^ = (the CP-conserving case) the A is a CP-odd Higgs boson {A' = A in 
such a case) and H, h are CP even. Such a statement is also true when ^ = tx j1 and/or 
As = Ae; as we shall see later in this section, for these particular values of parameters there 
is again no CP violation in the potential (P). 

For ^ = the Higgs boson masses can be calculated explicitly, and subsequently one can 
express the coupling constants Aj in terms of masses and a mixing angle defined through 

/cos« -smaW/^\ ^^Q^ 
n.2/ ysma cos ay V-"/ 

Let us now express the Ai, 2,3,4 in terms of the Higgs boson masses in the case ^ = (as we 
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have only four distinct masses, we leave the A5 as a free parameter). One gets 



A4 



A. 



As 
4 



Ai = -v~^ 
2 



22,22 ^aCjS j 2 2\ 



tan (i 



As 
4 



tan^ /? 



'Ill 



\ -'- -2 r 2 2 I 2 2 1. ni 2 2\1 As 

A2 = 2^ + - Sc.c/3tan/5(m^ - my^)\ - — 



tan^ /5 



Note also that the matrix of the quadratic form of the scalar fields is the Hessian of 
the potential at its minimum. The condition for the existence of a minimum is that the 
Hessian is positive definite, and this in turn means that the Higgs boson masses (squared) 
are positive. 

Finally, let us discuss briefly the particular cases ^ = 0,,^ = 7r/2 and As = Ae. The case 
^ = represents a model without CP violation within the scalar sector, as it is described 
in jS]. The case ^ = 71/2 can be analyzed easily in the parametrization ®; using this, the 
potential can be viewed as the case ^ = with the change of notation 



A, 



A, 



(12) 



Thus, the two cases are equivalent. When Ag = As, the ^-dependent part of the potential 
can be recast as 



As Re($l$2 



V1V2 



cos^ 



Ae ^Im($|$2 



V1V2 



sin^ 



Af 



(13) 



The remaining terms do not depend on the relative phase between $1 and $2? so that the 
phase factor e*^ can be transformed away and one thus again has a CP-conserving case. A 
particular consequence of such an analysis is that for z/ = there can be no CP violation. 



3 LQT method 

For finding the upper bounds on the Higgs boson masses we will employ the well-known 
LQT method invented three decades ago This method relies on imposing the condition 
of perturbative (in particular, tree-level) unitarity on an appropriate set of physical scat- 
tering processes. Within a renormalizable theory, the scattering amplitudes are "asymp- 
totically fiat", i.e. they do not exhibit any power-like growth in the high-energy limit. 
However, the dominant couplings are typically proportional to the scalar boson masses 
and one can thus obtain useful technical constraints on their values. In the pioneering 
paper |3] the method was applied to the minimal SM, and several groups of authors em- 
ployed it subsequently within models involving an extended Higgs sector, in particular the 
THDM (cf. IHj, [Zj, |H])- The results of various authors differ slightly, so it perhaps makes 
sense to reconsider the corresponding calculation and present, for the sake of clarity, some 
additional technical details of the whole procedure. 
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In the spirit of the LQT approach, our analysis is based on the condition of tree-level 
S'-matrix unitarity within the subspace of two-particle states. Instead of the unitarity 
condition used in the original paper jH], we can adopt an improved constraint for the 
s-wave partial amplitude A4o, namely 

iRe^ol < I (14) 

(cf. Pni)- Note that the tree- level matrix elements in question are real, and in the high- 
energy limit their leading contributions do not involve any angular dependence. Thus, the 
Aio generally coincides with the full tree- level (asymptotic) matrix element A^, up to a con- 
ventional normalization factor of IGir appearing in the standard partial-wave expansion. 
The effective unitarity constraint (fT^ then becomes 

|A^|<87r (15) 

For an optimal implementation of the unitarity constraints we will consider the eigenvalues 
of the matrix Mij = Aii^j where the indices i and j label symbolically all possible two- 
particle states. Having in mind our primary goal, we take into account only binary processes 
whose matrix elements involve the Higgs boson masses in the leading order, in particular 
in the 0{E^) terms. Invoking arguments analogous to those used in the original paper 
one can show that the relevant contributions descend from the interactions of Higgs scalars 
and longitudinal vector bosons. Using the equivalence theorem for longitudinal vector 
bosons and Goldstone bosons (see e.g. [Sj, JE]) one finds out, in accordance with the LQT 
treatment, that the only relevant contributions come from the amplitudes involving Higgs 
bosons and unphysical Goldstone bosons (that occur in an i?-gauge formulation of the 
theory). It means that we will examine the above-mentioned matrix Mij , including all 
two-particle states made of the scalars (both physical and unphysical) w^, z, H^, y4°, H, h. 
It is not difficult to see that the leading terms in the individual amplitudes are determined 
by the direct (contact) quartic scalar interactions, while the triple vertices enter second 
order Feynman graphs and their contributions are suppressed by the propagator effects in 
the high energy expansion. 

As noted above, we will be mainly concerned with the eigenvalues of the two-particle 
scattering matrix. It means that for our purpose we can consider, equivalently, any unitary 
transformation of the matrix Mij. In particular, it is more convenient to take, instead of the 
Mij, a matrix consisting of the scattering amplitudes between the two-particle states made 
of the "particles" w^,Za,ha corresponding to the parametrization (0). The eigenvalues of 
this matrix can be found in the earlier paper [H]. 

Matrix elements for the scattering processes corresponding to the two-particle states 
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{wiW2 , W2W1 , hiZ2, h2Zi, Z1Z2, hih2) form the submatrix 









"1^2 




Zl22 


^1^2 






/ 2Aq -I- ^ + ^ 


4 /As _ A6\ 


2A4 2"^^ 


2 ^4 2 6 


— A4 _|_ As 
2 2 


— A4 I 
2 


As \ 
2 \ 






2A3 + f + f 


^A4 + lAg 
4(f + ¥) 


IA4 — ^Ag 


— A4 _|_ As 
2 2 


— A4 _|_ 
2 


As 
2 


/ll22 


-Y^i + 1-^6 


IA4 — lAg 


As A^ 

4 (V + k) 










/l22l 


IA4 — ^Ag 


-2^A4 + lAg 


As As 
2 2 












— A4 _|_ As 
2 2 


— A4 _|_ As 
2 2 








4(f + ^) 

As Ag 
2 2 


As _ 
4(f + 


Ae 


hih2 


— A4 _|_ As 
\ 2 2 


— A4 _j_ As 
2 ~'~ 2 











with eigenvalues 



(16) 



(17) 



1 5 
ci = 2A3 — A4 — -A5 + -Ag 

62 = 2A3 + A4 - ^As + ^Ag 

5 1 
/+ = 2A3 — A4 + -A5 — -Ag 

1 1 
= 2A3 + A4 + -A5 — -Ag 

/i = /2 = 2A3 + ^As + ^Ag 

Another submatrix is defined by means of the states {wfwi, w^w^, ^^^); it 

reads 



Z1Z2 
hih\ 

^2^2 
V2 



4(Ai+A3) 
2X3 + ^ + ^ 

V2{Xi+X3) 

A 



V2{Xs + ^ 
\ V2{Xs + ^ 



2A3 + - 



and its eigenvalues are 



he 

2 



4{A2+A3) 
V2{Xs + ^) 
V2(A3 + ^) 
V2(A2+A3) 
V2{A2+A3) 



v^(Ai+A3) 
3(Ai+A3) 



2(^ 



2(^ + ^ 



V2{Ai+A3) 

V^(A3 + ^) 

A1+A3 
3{Ai+A3) 



V2 

V2(A3 + ^) 
y2(A2+A3) 

2(^ + ^) 
3(A2+A3) 



V2 



V2(A3 + ^) \ 
V^(A2+A3) 



2(^ + ^ 
2(^ + ^ 



2(^ 



(^ + ^) 2(^ + 4^) 2(^ + ^) 3(A2+A3) / 



(1^ 



a± = 3(Ai + A2 + 2A3) ± V9(Ai - A2)2 + [4A3 + A4 + \{\^ + Ag)]^ 



fei = Ai + A2 + 2A3 ± \l{Xi - X2f + K-2A4 + A5 + Ag 



c± = Ai + A2 + 2A3 ± y^(Ai - A2)2 + i(A5 - Ag)^ 
A third submatrix 

fti.i /2(A2 + A3) |(A5-Ag) 
h2Z2 \ 2(^5 — Ag) 2 (Ai + A3) 



(19) 



(20) 



has eigenvalues c± (see (HH)). Finally, there are submatrices corresponding to charged 
states {hiwf, h2wf, ziwf, Z2wf, /iiw^, /i2wj, Z1W2, 22102)'- 









21U1J 


h 


Z2Wl 








/2(Ai + A3) 


— A4 _|_ As 
2 2 

2 (A2 + A3) 







f A4 - 


f Ae 


\ 


h2W^ 


— A4 _|_ As 
2 2 


IA4 — 
2(Ai + 


f Ae 














^A4 + lAg 


A3) 


— A4 _l_ 
2 


As 

2 




Z2W^ 


\ -2^A4 + lAe 





— A4 1 
2 


As 
2 


2(A2 + 


A3) 


/ 



(21) 



hi to. 



h4«,+ /2(A3 + ^) 



h2wj 

1A4 I As 
2 2 



^ + ^ 2(A3 + ^^) 
":rA4 + -^Afi 



ZlW^ 



2A4 



f Ae 



22 Wo 



\ ^A4 + lAe 



2 (A3 + f ) 

-A4 
2 



2:2^2 
IA4 - f Ae \ 


— A4 I As 



(22) 



^ + ^ 2(A3 + ^)/ 

Their eigenvalues are the /_, 62, /i, c±, b± shown above and, in addition, 

Pi = 2(A3 + A4) - ^As - ^Ae 

Unitarity conditions (|T5|l for the eigenvalues listed above give the constraints 

|a±|, |6±|, |c±|, |/±|, |ei,2|, l/il, bi| < Stt 

Note that an independent derivation of these inequalities based on symmetries of the 
potential can be found in the papers [HI llOj. 



(23) 

(24) 
Higgs 



4 Independent inequalities 

However, the inequalities are not all independent. Indeed, it is not difficult to observe 
some simple relations as 

3/1 = pi + ei + /+ 

362 = 2pi + ei (25) 
3/_ = 2p^ + U 

and this means that the inequalities |/+|, |ei| < Stt imply |/i|, |e2|, |/-| < Stt. Further, 
the eigenvalues in the remaining inequalities can be rewritten as 

a± = 3Ai23 ± V^(3Ai2)2 + i(/+ + ei + 2pi)2 

b± = A123 ± Y/(Ai2)2 + ;^(/+ + ei-2pi)2 (26) 
c± = A123 ± v/(Ai2)2 + ^(/+-ei)2 
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where A123 = Ai + A2 + 2A3 and A12 = Ai — A2. In the case A123 > the inequahties for 
the a_, c_ follow from a^,b+,c^ < Stt. For A123 < the situation is similar, with 
interchanges (a, b, c)± (a, b, c)=p and A123 — A123. 

The authors [Zj noticed that among the latter inequalities, the strongest one is a+ < Sir; 
Indeed, using ()24|1 and one can show that for A123 > the remaining ones follow from 
it. In the case A123 < the same statement is true concerning a+ < Stt. 

Thus, it is sufficient to solve the inequalities 

|a±|,|/+|,|ei|,bi| <87r (27) 

In fact, the inequality a_ < Stt need not be taken into account in subsequent discussion; 
it turns out that this is weaker than the remaining ones and does not influence bounds 
in question (one can verify a posteriori that our solutions satisfy the constraints a_ < Sir 
automatically). 



5 Upper bounds for Ma and M± with ^ = 

Before starting our calculation, let us recall that the condition ^ = means that the Z2 
symmetry-breaking parameter u becomes u = (see dH}). To proceed, we shall first fix 
convenient notations. The LQT bound for the SM Higgs mass sets a natural scale for our 
estimates, so let us introduce it explicitly: 



""^^^ = ^1 w ^ ^ '^''^ ^^^^ 

(note that in writing eq. ()28|1 we do not stick strictly to the original value jS], using rather 
the improved bound ^Sj)- In the subsequent discussion we shall then work with the 
dimensionless ratios 

M = (29) 

"^LQT 

instead of the true scalar boson masses (denoted here generically as m). Further, an 
overall constant factor 167r/3 can be absorbed in a convenient redefinition of the coupling 
constants, by writing 

a: = ^ (30) 

Finally, we introduce new variables 

X = Ml + Ml Y = Ml-Ml Z=^-^Y (31) 

siri ^ jO 

that will help to streamline a bit the solution of the inequalities in question. 
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Using equations (fTTj) and the definitions shown above, the A' can be expressed as 



K 
K 

a; 



Y 



-AC 



Z 
1 



A' 



12 



A' 



123 



Ml 
Ml 

1 sin 2a, 
4 sin 2/3 

—1— [{X - 2\',) cos 2/3 - F cos 2a] 

2 sm 2/5 

1 A' 

5 — (X - r cos 2a cos 2/3 - 2 At) + ^ 

2sin2 2/?^ ^ ^^2 



(32) 



Let us now discuss the possible bounds for the M±, Ma- These can be obtained from the 
inequahties for |ei|, |/+|, \pi\, which read, in our new notation 



Mi 



^-K + 2Ml 



^Mi 



^Mi 



-M'i 



3 

< - 

- 2 

3 

< - 

- 2 

3 

< - 

- 2 



(33) 



The relations are linear with respect to the M^, and one can thus view the domain 
defined by these inequalities as a hexagon in the plane {M^, M^). Then it is clear that 
the highest possible value of a mass variable in question will correspond to a vertex (or 
a whole hexagon side). By examining all possible cases one finds easily that for M^, such 
a "critical" vertex satisfies the condition — /+ = pi = Svr; in view of (jHHj) this means that 
it corresponds to the values 



(M|,Ml) = (l + A'5,l + Z + 2A'5; 



(34) 



Such a maximum value of the is indeed formally admissible (in the sense that by 
reaching it one does not leave the parametric space of the considered model). To see this, 
one can substitute in eq. dSl = Ag, M^ = 1 + Ag, M| = 0, a = vr - /3. Thus, the 
bound becomes 

Ml<l + A'5 (35) 

Similarly, for the extremal solution corresponds to a hexagon vertex defined by ei = 
— /+ = Stt and its coordinates in the (M^, M^) plane are then 



:i + | + ^A'5,i + A'5) 



(36) 



(Ml, Ml) , 2 . 2- 

The parameter values that saturate this maximum are analogous and one has to take 
M| = A'5/2, M|r = 1 + A'5, M| = 0, a = vr - (3. In this way, the bound for Ml becomes the 
same as that for the M^ , namely 



Ml<l + X' 



(37) 
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6 Upper bounds for M/^, Mh with ^ = 



Let us now proceed to discuss the upper bounds for Mh and Mh. If we considered the 
relevant constraints without any further specification of the scalar bosons h and if, we 
would get the same result for both particles, since their interchange corresponds just to 
the replacement a —a (cf. eq. (fTUj) ). Thus, let us add the condition M^ < M^ 
(i.e. Y > 0). In such a case, we will solve just the inequality a+ < Sir (which puts the most 
stringent bounds on the variables X, Y) and in the obtained solution we will constrain the 
Ma, M± so as to satisfy the rest of the inequalities. 

The basic constraint a+ < Stt is quadratic with respect to the X, Y and reads (cf. the 
expression (HH))) 

(X-Fcos2acos2/?) - A5(2-sin2 2(3) + 

' \{X - 2A' ) cos 2(3 -Y cos 2a? + (-\ sin^ 2(3 ( _ ^ + A/2 + M)' < gin^ 2(3 

\3/ V sin2p 2 2 / 

(38) 

To work it out, we will employ the following trick: As a first step, we will consider a simpler 
inequality, which is obtained from (34) by discarding the second term under the square 
root; in other words, we will first assume that 

sin 2a At , Ml ^ 

Y ^ + Ml + — = (39) 

sin 2/3 2 =^ 2 ^ ' 

Of course, the "reduced" constraint 

X - r cos 2a cos 2p - \'^{2 - sin^ 2/5) + |X cos 2p-Y cos 2a - 2A5 cos 2/?| < sin^ 2/3 (40) 

is in general weaker than the original one. Nevertheless, in a next step we will be able to 
show that the obtained mass bound does get saturated for appropriate values of the other 
parameters (such that the condition ()39|1 is met) - i.e. that in this way we indeed get the 
desired minimum upper mass bound corresponding to the original constraint ()38p. Thus, 
let us examine the inequality (jlUj) . Obviously, we have to distinguish two possible cases: 

1. (X-2A'5)cos2/5 > Fcos2a. 
Then one has 

X(l + cos 2p) - Y{1 + cos 2/3) cos 2a - + cos 2/3) < sin^ 2p (41) 

Making use of our assumption, we can get from ()41|) a simple constraint that does 
not involve Y, namely 

X < 1 + A5 (42) 

(to arrive at the last relation, we had to divide by the factor 1 — cos 2/3; when it 
vanishes, we can use directly the original inequality PHj) and get the same result). 
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2. (X-2A'5)cos2/5 < Fcos2a. 

In a similar way as in the preceding case, the inequahty pUj) imphes the same bound 



Thus, having constrained X = Mfj + M| according to we can obviously also write 

Ml<l + A'5 (43) 

Now, it is not difficuh to see that for M| = 0, M| = A'5 + i, M| = 1 + A'5, a = vr - 
eq. is satisfied with M]j = 1 + A5 and means that PHj) represents the mass upper 
bound pertinent to the original unitarity constraint ()38|). 

The bound for the Mh is obtained from ()42|1 by using there our subsidiary condition 
Mh < Mh] one thus has 

Ml<lil + K) (44) 

The upper limit in ^ gets saturated (i.e. M| = |(1 + A5)) for Mh = Mh, Ml = 0, 
M^ = A5/2, a = 37r/4,/3 = 7r/4. It is worth noticing that here we have fixed a particular 
value of the angle f5 , while all previous constraints were independent of (3 (i.e. for any j3 
we were then able to find an appropriate value of a). A more detailed analysis shows that, 
in general, the upper bound for the M^ indeed depends explicitly on the (3. To derive the 
corresponding formula, we consider the boundary value Mh = Mh (i.e. Y = 0) and use 
also eq. (jHSI)- The inequality (PH|1 then becomes 

Ml - K (1 - ^) + |M^os2/3 - A',cos2/3| < ^ (45) 

To work it out, we will assume that M| > A5 (taking into account pij) this means A5 < 1; 
in fact, one can do even without such a restriction, but for our perturbative treatment only 
sufficiently small values of the A5 are of real interest). The inequality then becomes 

m2 / (l-A-5)(l + cos2/3)(l-cos2/3) 
^'^ - ~~2 l + |cos2/3| + 

Obviously, the maximum bound ()44|1 is recovered from the last expression for (3 = 71 /A. 
Let us also remark that the choice a = ir — j3 comes, as in all previous cases, from the 
requirement Z = —Y. 



7 Upper bound for the lightest scalar for ^ = 

One can notice that any Higgs mass upper limit discussed so far gets saturated only 
when at least one of the other scalar masses vanishes. Thus, another meaningful question 
arising in this connection is what can be an upper bound for the lightest Higgs boson 
(within a considered set of the five scalars h, H, A^, H^). Let us first take h to be the 
lightest scalar state; it means that in our analysis we will include the additional assumption 
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Figure 1: The region of admissible values of if h is assumed to be the lightest scalar 

Mh < Mhi Ma, M±. The procedure we are going to employ is a modest generalization of 
the earher calculation [?]. Squaring the inequality (|HHj) one gets 

(X - X,f sin^ {Y - Y,f > (47) 

where Xq,Yq and R depend on X'^,a, P, + M\/2. This inequality defines the domain 
bounded by the hyperbola shown in Fig. but the original constraint (jHH|) corresponds 
just to its left-hand part. In order to find the solution, one should realize that the slope of 
the asymptote with respect to the X-axis must be greater than the slope of the straight 
lines X = ±Y (this follows from the fact that the coefficient 1 — | sin^ 2a, multiplying the 
in (j37j), is less than one). Because of that, the maximum value of the Mh corresponds 
to y = and a+ = Sir, and we are thus led to the equation 



X - X'rp - sin^ 2/5) + J cos2 2(3{X - 2X'^y + - sin^ 2/3 ( M| + ^ - ^ j = sin^ 2(3 (48) 

It is clear that for smaller M±, one has a bigger value of the Mh, so the needed 
upper estimate is obtained for M± = Ma = Mh (note also that from Y = one has 
X = 2M^). In this way one gets an equation for maximum Mh- 



2Ml - \'^{2 - sin^ 2(3) + ^ A{Ml - X'^Y cos^ 2(3 + sin^ 2(3{Ml - -X'^f = sin^ 2/3 (49) 

From eq. ()49|1 one can calculate the M| as a function of sin^ 2/3. It can be shown that for 
Ag < 3/5 this function is increasing, i.e. the maximum is reached for (3 = n/A and its value 
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Figure 2: Dependence of lightest boson mass on (3 



becomes 



Ml 



(50) 



We do not display the explicit dependence of the maximum on the but it is clear that 
the solution of eq. ()49|) is straightforward. Finally, we should also examine the cases where 
the lightest Higgs boson mass is either or M± . However, from the above discussion it 
is clear that both these extremes occur when Mh = Ma = M±. 

Similarly, from eq. ()48|1 one can derive a constraint for the mass of the lightest neutral 
scalar boson (which we denote M„). In this case we substitute there X = 2M^, M\ = M^, 
= and obtain thus the equation 



2Ml - X'^{2 - sin^ 2(3) + \ U{M^ - X'^y cos^ 2(3 + sin^ 2/3 



A/2 



sm^2p (51) 



From eq. (jKT|l one then obtains the as a function of sin^ 2/?, which is increasing for 
A5 < 1. Its maximum reached at /3 = 7r/4 becomes 



(52) 



8 Numerical solution for ^ 7^ 

In the general case with 7^ (i.e. with CP violation in the scalar sector) we have not 
been able to solve the inequalities (jTTj) analytically, so we had to resort to an appropriate 
numerical procedure. The main result we have obtained in this way is that for small values 
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of the parameter v (see eq. (0)), in particular for v' G (0, 0.3) , the upper mass bounds in 
question are the same as for ^ = 0. The interval has been chosen such that the variations 
in the upper estimates be at the level of 50 — 100%, the validity of our theoretical estimates 
is guaranteed up to i^' < 3/5 (see the remark below eq. P^jl ). 

Our numerical procedure consists in solving the inequalities fl27|l on the space of pa- 
rameters A'i 2 3 4 5 6 ^^'^ ^ restricted by the condition Q, where one also adds constraints 
for the existence of a minimum of the potential (^: A4 > (i.e. > 0, see flllj) ) and the 
requirement of positive definiteness of the matrix Q (i.e. fn^^i^ > 0). On this parametric 
subspace we have looked for the maximum values of the following quantities: 

1. Mass of the charged Higgs boson m± (see Fig. E)) 

2. Mass of the lightest Higgs boson (see Fig. Hj) 

3. Mass of the lightest neutral Higgs , i.e. the lightest one among the A, H, h (see Fig.^)) 

4. Mass of the heaviest neutral Higgs, i.e. the heaviest among the A, H, h (see Fig. IH)). 

Let us remark that in this case we have not distinguished between A and h, H, which are 
superpositions of the CP-odd and CP-even states. 

In our plots we display, apart from the dependence of masses in question on the u, also 
the values of the parameter ^ in the case A5 = Ag and A5 7^ Ag respectively, in order to be 
able to distinguish the extreme cases without CP violation (^ = /c7r/2 or A5 = Ag, see the 
discussion in Section |21). From Figs. EHHEliniit can be seen that all examined mass upper 
bounds are reached just in the aforementioned extreme cases. In view of this, we can make 
use of our previous analytic expressions, except for the case El which we have not solved 
analytically. 

Our results have been simulated by means of the computer program Matlab 6.0, package 
optim, with the help of the function fmincon. The numerical errors are mostly due to an 
insufficiently smooth condition for the positive definiteness of the matrix 0. 

9 Conclusions 

In the present paper we have reconsidered upper bounds for the scalar boson masses within 
THDM, by using the well-known technical constraint of tree-level unitarity. Our analysis 
should extend and generalize the results of some previous treatments, in particular those 
obtained in the papers [Hj and Although we basically employ the traditional methods, 
we have tried to present some details of the calculations not shown in the earlier papers — 
we have done so not only for the reader's convenience, but also to provide a better insight 
into the origin of the numerical results displayed here. As we have already noted in the 
Introduction, some new relevant papers on the subject have appeared quite recently (see 
jni Cni IH] ) • in these works, the structure of the unitarity constraints is discussed in detail 
within a rather general THDM, but there is no substantial overlap with our results, since 
our main point is rather a detailed explicit solution of the inequalities in question. 
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Figure 3: Charged Higgs boson, with theoretical estimate (jH^ 
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So, let us now summarize briefly our main results. We have found upper limits for 
Higgs boson masses in dependence on the parameter v that embodies an information 
about possible flavour- changing neutral scalar-mediated interactions. The upper bounds 
are seen to grow with increasing v (see Tab^. On the other hand, this parameter cannot 
take on large values (to avoid a conflict with current phenomenology), and thus it makes 
no real sense to consider the mass estimates for an arbitrary v ; in the present paper we 
restrict ourselves to < 0.4 (cf. the condition used when deriving the relation (jSU))). In 
the case with no CP violation in the scalar sector (,^ = 0), the relevant results are obtained 
from the inequalities (jHHjl . (jSZ|), PHjl - and the bound for the lightest scalar is shown 
in eq. ljHnjl (where one should also pass from A5 to A5 according to (jHUjl ). In Section |H1 we 
have then verified that in the CP-violating case these values remain the same. The results 
are shown in Tab. Q where we have singled out the case 1/ = that corresponds to the 
absence of flavour- changing scalar currents. Let us remark that in the CP-violating case 
we do not distinguish between the H and A, and in the CP-conserving case the bounds 
for H and A are the same. 

Further, we have calculated an explicit dependence of the upper limit for the M/^ on 
the angle /? in the case with ^ = 0. The analytic expression reads 



(cf. fl4fij) with the A5 retrieved). The dependence of the relevant bound for a lightest scalar 
boson can be obtained from eq. (imj) and the results for some particular values of the A5 
are depicted in FigO 

For = and ^ = 0, our results can be compared directly with those published in [Zj. 
We get somewhat stronger bounds for iua and m± since, in addition to the set of constraints 
utilized in we have employed also the inequality pi < Svr, which stems from charged 
processes (cf. the end of Section not considered in [7]. On the other hand, our estimates 
for m u , vTLh and the lightest scalar coincide with the results [7] , since the above-mentioned 
extra inequahty is not used here. It is also noteworthy that the upper hmits for and 
mn coincide with the SM LQT bound if they are estimated separately and, depending 
on the number of the simultaneously estimated Higgs scalars, the coefficient 1/2 appears 
when we take two of them and 1/3 when all of them are considered. 

In the case ^ = and A5 7^ comparison with [S] is possible. Here we can compare 
only the corresponding numerical values, which turn out to be approximately equal when 
A5 = 0. However, for A5 = our results obviously differ from those of jH]: in particular, the 
bounds for mA, Tn± displayed in [Hj appear to decrease with increasing A5. The authors [Hj 
state that they used some fixed values of the angle /5; for the purpose of a better comparison 
we have therefore calculated the /3-dependence of the upper bound for m/^, with the result 
shown in As it turns out, the and m± do not depend on (5 in this case. 

Finally, let us mention that in the CP-violating case we have not been able to get 
analytic results; we have only shown, numerically, that the maximum values of the masses 
in question are obtained for ^ = 0, i.e. the upper mass bounds are the same as in the case 
with no CP violation in the scalar sector. 




(53) 
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